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Abstract 

The level splitting formulae much discussed in the study of spin tunneling in 
macroscopic ferromagnetic particles and previously derived only by compli- 
cated pseudoparticle methods for the ground state, are derived from those of 
eigenvalues of periodic equations and extended to excited states. 
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In recent years the tunneling of quantum spins and the possibihty of 
its observation in macroscopic ferro- or antiferromagnetic particles has at- 
tracted considerable interest as the rapidly expanding literature in the field 
demonstrates [|l|, 0. The principal idea is that in such particles of mesoscopic 
size the electronic spins can form an aligned magnetic state which can as- 
sume several directions so that quantum mechanics suggests the possibility 
to lift this degeneracy by tunneling from one direction to another. This tun- 
neling effect and its possible suppression in half-integer-spin particles (with 
or without applied magnetic fields) has been explored in numerous investiga- 
tions 1^, H, 1^ . In most of these the first step is the conversion of the discrete 
spin system into a continuous one |^ by transforming spin operators to 
canonical operators p, q satisfying the quantisation relation 

[Q.p] = -s 

where s is the spin quantum number of the particle. The calculation of the 
basic semiclassical Hamiltonian is quite involved and can be found in the 
literature ^, and will therefore not be repeated here. Once the approxi- 
mate effective quantum hamiltonian has been constructed, the next step is 
ususally to proceed to the appropriate Lagrangian and then to the applica- 
tion of the path-integral method. Here, however, a problem arises. Since the 
effective Hamiltonian is that of a particle with position- or field-dependent 
mass (see below) , the path-integral method must start from the phase-space 



path-integral, as pointed out long ago [jTO|- The resulting evaluation of the 
path-integral by expansion about a classical configuration P, [ll|, |l^ which 
itself is already a complicated expression in many cases (involving sometimes 
elliptic functions) is so complicated that one would like to have a verification 
of the resulting formulae, at least in some limiting cases, although, no doubt, 
the pseudoparticle methods applied to these theories are of considerable in- 
terest in themselves. Nonetheless, particularly experimentalists interested 
in spin tunneling phenomena might welcome a more straightforward deriva- 
tion of the relevant level splitting formulae. Our objective here is therefore 
to point out that in certain cases such a simpler verification is possible by 
reference to some standard (though possibly not so well-known) quantum 
mechanical problems, and yields even the splitting of the appropriate ex- 
cited state levels (a definite advantage over the instanton method). In the 
following we therefore demonstrate that in the case without an applied mag- 
netic field the level splitting formulae for the ground state as well as excited 
states can be obtained from those of Schrodinger equations with periodic 
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boundary conditions. In the limit of vanishing position- or field-dependence 
of the effective mass the comparison equation is the Mathieu equation, and 
in the case with position-dependence it is the (even less well-known) Lame 
equation. In both cases the results are easily extendable to excited states 
and the range of validity of the results can be clearly inferred from the valid- 
ity of the asymptotic solutions of these equations. The results, which verify 
previous path-integral calculations, also demonstrate that the field depen- 
dence of the effective mass provides only a minor change. We expect that 
the derivation given here may be helpful also in the verification of the results 
of some other tunneling model theories. 

The Hamiltonian of the theory describing the ferromagnetic particle is 
given by 

H = K^Sl - K2SI 

in terms of spin operators Sx,Sz and Ki > K2 > 0, the two contributions 
describing XOY-easj plane anisotropy with easy axis along the x direction 
[§]. As mentioned above the spin operators are transformed to canonical 
operators p, </> with Sz = sp and = -. One then obtains an effective 
Hamiltonian which has been used in the following form in the literature 

B 1, 



where 



and —IT < < vr. We consider first the case of constant mass, i.e. A = 0, 
although in actual fact [§, || A = K2/K1. The problem is to determine 
the splitting of asymptotically degenerate energy levels of the Schrodinger 
equation Hip = Eip for the given periodic potential. Setting 

the Schrodinger equation becomes the Mathieu equation 

ilj'^cj)) + (a - 2h^ cos 20) V = (4) 

The tunneling splitting of the eigenvalues A of this equation (which is twice 
the shift of a single level) has been derived in the literature [|13| from periodic 
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boundary conditions. In terms of the quantum number go = 2n + 1, n = 
0, 1, 2, ... the dominant term can be written 

" (8,r)4|i(*-l)|! ^ 'h'J 

For the ground state (i.e. go = 1) this yields the level splitting of our original 
problem here, i.e. 




(6) 

Comparing this with the result of p| in the same limit (i.e. for no applied 
field and the parameter a — 0) in [Q) one obtains the same result except 
that in |^ the exponential factor is missing. The exponential factor is crucial 
evidence of the nonperturbative nature of the calculation (in the calculations 
underlying (5) the effect of boundary conditions), and its argument usually 
represents the euclidean action of the appropriate instanton. The exponen- 
tial factor is missed in 0, presumably as a consequence of the approximation 
considered which then is the leading approximation of that under discussion 
and is therefore indicative of the approximation there employed. This also 
shows the conditions on the parameters under which the result is valid. The 
parameter h"^ of eq.(|^) has to be large, and the result is the dominant con- 
tribution of an asymptotic expansion in descending powers of h"^. This con- 
dition implies that >> 46, as also realised in |^. An additional advantage 
of the Schrodinger equation method given here is that it yields immediately 
also the level splitting for excited states, which cannot be obtained with the 
usual instanton method, but instead requires the corresponding consideration 
of nonvacuum (or periodic) instantons or some other method as in ||12|| . 

We now want to take the 0-dependence of the effective mass m(0) into 
account. The expression itself in eq.(|^) is reminiscent of an elliptic expression. 
This suggests searching for an analogy with an elliptic equation like the Lame 
equation. We set x = sin0 and = {1 — a;^)~4$. We then expand the 
coefficient of $ in the resulting equation in ascending powers of x and obtain 
an equation which we can write (with — = — + Ax^) 

(1 - x^){l - Ax2)$" + |1 [2 + 2A + x\3 - 6A) + x\3 - 3A) 
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+ 



\E-s^K^) y 

Ki 2, 



1 Ax' + 0(x^A2) $ = (7) 



The Lame equation is best known in the form WA 



w" + |A - n{n + l)k^sn^z^ w = 



(8) 



The eigenvalues A and eigenf unctions w of the equation with periodic bound- 



ary conditions can be found in the hterature [0]. The level splitting analo- 
gous to that of the Mathieu equation is not so well-known but also available 



in the literature |T5|. In fact, under certain conditions the Lame equation 
reduces to the Mathieu equation. In eq.(^ n is an integer in the case of 
Lame polynomials, and k is the elliptic modulus of the Jacobian elliptic 
functions. The quantity A is the eigenvalue. To a good approximation we 
can convert eq. @ into one similar to eq. (J^) . To this end we set t = snz and 

w = [{1 — t^)(l — kH"^)]"^ S. The equation for S is then 







(9) 



Comparing eq.(|^ with eq.(|D we can identify the following quantities 

A (10) 



The level splitting of the eigenvalue A of the Lame equation is known to 
be in leading order of = n{n + l)k'^ |]T5[ 



AAgo = Ak 



2\2 fl+k' 
Vl-A; 



TT 



_8k\ 2 



5 go 



l + 0(- 



(11) 



where for small values of k 

l + k 



In 
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In the case of the ground state (go = 1) eg. ([I0|) imphes the sphtting of 
the Schrodinger eigenvalue E given by 




(12) 



We compare this result with that given in [§] (there eq. (9a)) and calculated 
with the help of a complicated instanton method. For the comparison we 
have to set in our result k'^ = ^, so that eq.(0) becomes 

16Ki 



AEi 



7S V^Ti^' ^ 



(13) 



The corresponding expression obtained in the literature M by the instanton 
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1 + 4)e ^'^ 
46 



(14) 



We observe that except for small differences in the coefficients of the effects 
of the field dependence of the effective mass, the expressions obtained by 
these totally different methods agree with one another. We interpret the 
small differences in the coefficients as due to unavoidable and somewhat 
different approximations in the two methods. Our procedure here extends the 
previous results in several ways. First, one can clearly see that the result is 
the dominant contribution of an asymptotic expansion in descending powers 



of K 



This is therefore the parameter which is assumed to be large, in 



agreement with the condition on h"^ above and experimentally relevant values. 
Secondly, apart from the simple reference to the result in the literature, this 
Schrodinger method yields also immediately the level splitting of excited 
states which are much more difficult to obtain with a pseudoparticle method. 
Since these are easily obtained from the cited formulae above we do not write 
them down in detail for the problem under discussion. 

Concluding we can say that the level splitting formulae, which have been 
much discussed in connection with spin tunneling in macroscopic particles 
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but have previously been derived only by complicated pseiidoparticle path- 
integral methods can readily be obtained in equivalently good approximation 
from the level splitting of eigenvalues of periodic equations, which admittedly 
are not so well known but have been studied extensively in the literature. 
This simpler way of deriving these formulae should appeal particularly to 
experimentalists. 
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